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The aerodynamic performance of a high-load low-pressure turbine blade cascade has been analyzed for three dif- 
ferent distributed surface roughness levels (Ra) for steady and unsteady inflows. Results from CFD simulations 
and experiments are presented for two different Reynolds numbers (300000 and 70000 representative of take-off 
and cruise conditions, respectively) in order to evaluate the roughness effects for two typical operating conditions. 

Computational fluid dynamics has been used to support and interpret experimental results, analyzing in detail 
the flow field on the blade surface and evaluating the non-dimensional local roughness parameters, further con- 
tributing to understand how and where roughness have some influence on the aerodynamic performance of the 
blade. The total pressure distributions in the wake region have been measured by means of a five-hole miniatu- 
rized pressure probe for the different flow conditions, allowing the evaluation of profile losses and of their de- 
pendence on the surface finish, as well as a direct comparison with the simulations. 

Results reported in the paper clearly highlight that only at the highest Reynolds number tested (Re=300000) 
surface roughness have some influence on the blade performance, both for steady and unsteady incoming flows. 
In this flow condition profile losses grow as the surface roughness increases, while no appreciable variations have 
been found at the lowest Reynolds number. The boundary layer evolution and the wake structure have shown that 


this trend is due to a thickening of the suction side boundary layer associated to an anticipation of transition 


process. On the other side, no effects have been observed on the pressure side boundary layer. 
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Introduction 


The first experimental work on the effects of rough- 
ness was carried out by Nikuradse [1], who studied the 
losses in sand-roughened pipes varying the Reynolds 
number and the roughness itself. In order to characterize 
the roughness height, he defined a non-dimensional pa- 
rameter (k, ) that compared the roughness with the visc- 
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ous length scale of the boundary layer. He found that, in 
the case of fully turbulent flows, the roughness didn't 
have any effects on losses for k,' < 5: in this case rough- 
ness peaks are entirely confined within the viscous sub- 
layer and the oscillations generated by the roughness are 
damped by the viscosity of the fluid. Nikuradse called 
this regime "hydraulically smooth". When 5 < k,* < 70 
both Reynolds and roughness had an influence on losses 
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Nomenclature 
C chord XYZ cascade axial, tangential and spanwise 
coordinates 
c f - U mean veloci 
” static pressure coefficient = Pie. y 
Pr P2 
c . = Ubar bar rotational speed 
4 total pressure coefficient = 242+ sar k 
Pn P2 
ks equivalent non-sand roughness ur friction velocity 
ks non-dimensional roughness parameter Greek letters 
C axial chord a flow angle 
d bar diameter v kinematic viscosity 
fia bar passing frequency o flow coefficient= U Upar 
f reduced frequency=fjq-C/U o total pressure loss coefficient 
G pitch Op profile total pressure loss coeffi- 
n Pn Pt2,m 
cient = ———— 
Pu P2 
h blade height Üa wall shear stress 
p static pressure Subscripts and superscripts 
Pi total pressure ~ pitch-averaged quantity 
1 at the cascade inlet 
2 at the cascade outlet 
engine 
a mass averaged exit total pressure = ref reference condition (midspan, steady inflow, 
Pro 
h/2g h/2g Re=300000) 
f fuxpadyaz/ f fuxdyaz 
0 0 0 0 
Pim mass averaged exit total pressure at midspan= x in the cascade axial direction 
? g g 
fuxpady/ fay 
0 0 
Ra distributed surface roughness w Wall 
Re Reynolds number =U;C/v 


because viscosity is no longer able to damp out the tur- 
bulent eddies created by the surface asperities (Flack and 
Schulz [2]): for this reason, the regime was considered as 
"transitionally rough". At last, for k, > 70 the losses only 
depend on the roughness level: roughness elements gen- 
erate turbulent eddies regardless of Reynolds number. In 
this condition the higher the roughness the higher the 
losses. This is the "fully rough" regime. 

On the other side, in the case of laminar pipe flows 
roughness has not been found to influence the pressure 
losses, as long as roughness is not able to trigger transi- 
tion. The possibility that roughness triggers transition in 
laminar flows depends on the roughness height compared 
with the boundary layer thickness. An accepted criterion 


to estimate the occurrence of roughness-induced transi- 
tion is that proposed by Braslow [3], who suggested the 
threshold of k,” = 19. 

However, the effects of roughness in real applications, 
such as over turbine blades, are further complicated by 
pressure gradients, free-stream turbulence intensity and 
passing wake unsteadiness that strongly affect the transi- 
tion process on both the pressure and the suction sides of 
the blade (e.g. Hodson and Howell [4]). Moreover, since 
high pressure and low pressure turbine blades are cha- 
racterized by sensible different loading distributions as 
well as operating Reynolds numbers, roughness effects 
may sensibly differ for these two modules. For low pres- 
sure turbines (LPT), operating at low Reynolds number 
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during cruise conditions, the effect of “as-cast” rough 
blades can reduce or eliminate separation bubbles via the 
"roughness-induced transition", thus reducing losses (e.g. 
Hodson and Howell [4], Vazquez and Torre [5]). Vera et 
al. [6] noted that the optimum surface finish of LPT 
blades does not necessarily correspond to the aerody- 
namically smooth surface, and that it can be possible to 
reduce the manufacturing cost removing the need to 
polish the blades. On the other hand, for attached boun- 
dary layers developing at high Reynolds number (such as 
in high-pressure turbine blades) the effect of roughness 
can provoke a significant increase of both losses and heat 


transfer coefficient as reported by Boyle and Senyitko [7]. 


Their tests show that roughness nearly doubled the loss 
level at high Reynolds number. 

From the numerical point of view many corrections 
have been proposed to account for wall roughness with 
the most used turbulence models. Aupoix and Spalart [8] 
formulated some extensions of the Spalart-Allmaras tur- 
bulence model to account for non-hydraulically smooth 
surfaces. More recently Aupoix [9] critically reviewed 
the corrections available for the k-w shear stress transport 
(SST) model. Other works can be found in literature de- 
scribing the combined effects of roughness with other 
dimensionless parameters such as the free-stream turbu- 
lence level, as documented in the LES computations of 
Nagabhushana Rao et al. [10]. Bellucci et al. [11] devel- 
oped a CFD roughness model and tested it against expe- 
rimental data relative to HPT steam turbines. 

An excellent review of the effects of roughness on gas 
turbines is reported in Bons [12]. 

However, the real operation of both HP and LP turbine 
blades is sensibly affected by the rotor-stator aerody- 


namic interaction (see Hodson and Howell [4], Satta et al. 


[13]), since wakes shed by the upstream row alter signif- 
icantly the boundary layer transition process, thus the 
influence of roughness. However published investiga- 
tions of the roughness effects in unsteady cases are quite 
few (e.g. Montomoli et al. [14]). 

The present activity is part of a joint research effort 
between GE Avio, the University of Genova, and the 
University of Florence, aimed at developing optimum 
design criteria for high-lift low-pressure turbine blade 
cascades, through the detailed understanding of the 
physical phenomena interesting this kind of blades. 

In this work the effects induced on the aerodynamic 
performance of a high-load LPT blade cascade by three 
different levels of surface roughness have been analyzed. 
Experiments have been carried out for steady and un- 
steady inflows at low and high Reynolds number: Re = 
70000 and Re = 300000 representative of cruise and take- 
off conditions, respectively. Blade aerodynamic loadings 
and total pressure distributions in the wake region have 
been measured for each roughness level, allowing the 
evaluation of profile losses as a function of surface rou- 
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ghness for the different flow conditions. Computational 
fluid dynamics has been used to analyze in detail the en- 
tire flow field on the blade surface for steady state condi- 
tions. CFD data also allowed to evaluate the distributions 
of the non-dimensional local roughness parameter k,". 
This helped to understand how and where roughness 
have some influence on the boundary layer evolution, 
thus further contributing to the understanding of the 
roughness effects on the aerodynamics performance of 
low pressure turbines. 


Experimental Facility and measuring techniques 


The experimental investigations have been carried out 
in the blow-down wind tunnel installed at the Aerody- 
namics and Turbomachinery Laboratory of the University 
of Genova. The test section is constituted of a 7 blade 
large scale planar cascade, representative of LPT blade 
profiles. The airfoil has been designed by GE Avio and 
scaled for low speed testing following the procedure de- 
scribed in [15]. Blades are characterized by a chord of 
120 mm and an aspect ratio of AR = 2.5 to ensure 
two-dimensional flow at midspan. The main geometrical 
parameters of the tested profiles are summarized in Table 
1. The central blade of the cascade has been tested in 
three different surface roughness configurations: “as 
casting”, “as drawing” and “super-polished”. These are 
representative of a very row, a standard and a very accu- 
rate degree of surface finish attainable with different 
manufacturing processes in the real engines. Fig. 1 shows 
a scheme of the test section. 


SLOT FOR BARS 
MOVEMENT 


DOWNSTREAM 
MEASURING PLANE 


FRONTAL BLEEDS 
TAILBOARDS 


Fig.1 Sketch of the test section 


Table 1 Cascade main geometrical parameters 


Chord length C=120mm 
Aspect ratio h/C = 2.5 
Pitch to chord ratio g/C = 0.685 
Inlet metal angle at; = 43.95° 
Outlet metal angle a = -65.8° 
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A turbulence generation grid has been installed up- 
stream of the cascade in order to generate an inlet turbu- 
lence intensity of 5.2%, which is closer to the real engine 
value. 

For what concerns the unsteady tests, a moving bars 
system has been adopted to simulate upstream wakes. It 
consists of a tangential wheel of radial rods rotating in a 
plane parallel to the cascade leading edges plane at a dis- 
tance of 33% of the blade chord upstream of it. The flow 
coefficient and the reduced frequency for the unsteady 
inflow tests were chosen to be representative of real en- 
gine operative conditions (y = 0.675 and f° = 0.69). The 
bar diameter was chosen so that the wakes shed from the 
bars produce the same losses as those generated by a 
typical upstream LPT row. Further details on the bars 
system are provided in Satta et al. [16]. The flow deflec- 
tion through the bars system due to the velocity defect 
associated with incoming wakes has been compensated 
by means of a counter rotation of the cascade. This action 
ensures the same mean incidence angle for the steady and 
the unsteady tests, as shown in Simoni et al. [17]. 

To obtain the different values of surface roughness, the 
central blade of the cascade was sandblasted with differ- 
ent sand sizes. The cascade aerodynamic performance 
has been measured for each surface roughness configura- 
tion for different Reynolds numbers (Re = 70000 and Re 
= 300000), under both steady and unsteady inflows. 

The blade aerodynamic loading has been surveyed on 
the two blades closer to the central one in order to ensure 
the identity of the boundary conditions during the analy- 
sis of the three different roughness configurations and to 
set the periodicity condition by means of adjustable tail- 
boards and sidewall inlet bleeds. To this aim these two 
blades were instrumented at midspan with a total of 22 
pressure taps connected to a Scanivalve. 

The time-mean total pressure distributions at the cas- 
cade exit plane have been measured by means of a 
five-hole miniaturized pressure probe, with 3 mm di- 
ameter head. The small probe diameter jointly with the 
large scale dimension of the cascade allows obtaining 
high accuracy measurements in the wake region. Mea- 
surements have been performed at midspan in a tangen- 
tial plane located at an axial distance of 30% of C, from 
the cascade trailing edge plane. 26 equally spaced points 
have been adopted to survey 1 pitch. A three-axis com- 
puter controlled traversing mechanism with a minimum 
linear translation step of 8 um has been employed to al- 
low high movement precision and spatial resolution. 
Pressures were measured by means of high-sensitivity 
high accuracy low-range (+ 620 Pa) SETRA differential 
transducers, ensuring accuracy better than +0.075 % of 
the transducer full-scale range. For each measuring point, 
10000 samples have been collected with a sampling fre- 
quency of 1 kHz. In order to reduce experimental errors 
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due to the inherent uncertainty of the pressure transduc- 
ers or to the change of environmental conditions, the 
time-mean total pressure flux leaving the cascade has 
been measured three times, and results averaged. The 
uncertainty in the evaluation of the total pressure loss 
coefficients is lower than 0.02 %. 


Computational Framework 


The TRAF code (Arnone [18]) was selected for the 
numerical investigations. The code solves the unsteady, 
three-dimensional, Reynolds-averaged Navier-Stokes 
equations in the finite volume formulation on multi-block 
structured grids. Convective fluxes are discretized by a 
2™ order TVD-MUSCL strategy build on the Roe’s up- 
wind scheme. A central difference scheme is used for the 
viscous fluxes. The concept of artificial compressibility 
of Chorin [19] is used to handle incompressible fluids by 
a time-marching approach. 

Wilcox’s low-Reynolds number k—œ model [20] is ap- 
plied in combination with the laminar kinetic energy 
(LKE) model, which enables to take into account the 
pre-transitional rise of the fluctuating kinetic energy 
(Mayle and Schultz [21]). The surface roughness model 
of Wilcox [22] was used, at the wall the turbulent kinetic 
energy is set to zero and the specific dissipation is mod- 
ified in order to enhance turbulence compared to a 
smooth wall. A detailed description of the three-equation 
turbulence/transition model and the correlation used to 
sensitize the LKE transition model to roughness can be 
found in Pacciani et al. [23] and Bellucci et al. [11]. The 
computational domain was discretized using an O-H grid 
for a total of about 70,000 cells. The mesh was selected 
on the basis of previous grid-dependence analyses [11]. 


Reproduction of roughness level for cascade 
testing 


The surface roughness of a turbine blade directly de- 
pends on the surface finish obtained with the manufac- 
turing process as well as on the deterioration to which the 
blade surfaces are subjected during their operation. The 
scaling of the roughness values from the real engine to 
the test case was performed taking into account the Ni- 
kuradse’s roughness non dimensional parameter k. In- 
deed, according with the technical literature this is the 
most suitable parameter to describe the effects of rough- 
ness on the boundary layer development. From Niku- 
radse’s definition of k,”: 

kt = kir 
v 

where, according to Schlichting [24], k, describes an 
equivalent "non-sand" roughness that has the same ef- 
fects on losses of sand-grains with diameter k,. With re- 
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gard to gas turbines, there are several correlations to 
evaluate k, [12], and it is difficult to determine the equiv- 
alent sand grain roughness for a given surface since the 
result may vary by up to a factor of five. Most of them 
give a liner relationship between k, and the measured 
values of distributed surface roughness (Ra). For the 
present case of sandblasted surface, probably the most 
appropriate proportionality coefficient is 8,9 proposed by 
Schaffler [25] or 5,863 proposed by Adams et al. [26]. 
Multiplying and dividing for characteristic velocities and 
length scale, and expanding the friction velocity it is 
possible to write k,” as: 


pt aks Ure Ue 
* C v U;,\ pU; 


Finally, recognizing some fundamental non dimen- 
sional parameters of the problem, it is possible to write: 


Akre (£ Re, Le, sc] 


Therefore k,' a w the effects of e on the 
boundary layer, directly depends on the cascade outlet 
Reynolds number, on the profile velocity distribution, the 
skin friction coefficient and on the ratio between k, and 
the blade chord. Since the friction coefficient can be con- 
sidered mainly a function of the Reynolds number and 
the loading distribution, the previous relation can be fur- 
ther simplified to: 


k U 
kt = f-| ,Re,,—* 
f r| pae 


This makes evident that the same dimensionless ratio 


k= 


S 


Cp = 


€ should be fixed, once working with blade loading 


and Reynolds number similarities, between the test and 
the engine conditions. It also highlights that the higher 
the Reynolds number the higher k,*. Thus, a surface of 
given roughness could be hydraulically smooth at low 
Reynolds number while at higher Reynolds roughness 
could significantly affect losses, as revealed by results 
shown in [7]. In order to achieve a proper assessment of 
roughness effects on the performance of LPT blades, 
three different surface roughness levels have been tested. 
The values of the surface roughness were chosen to be 
representative of those that can be achieved with differ- 
ent LPT blade manufacturing processes (expressed in 
terms of the root mean square of the peak/valley height, 
namely the Ra level): 

@ "super polished" blade (Rajesi; = 0.95 um), re- 
presenting a higher degree of finish as compared 
with the actual one (Raengine= 0.3 um), very expen- 
sive. 

"as drawing" blade (Rajesi; = 5.7 um), representing 
p degree of finish normally requested for the 
manufacturing of a LPT blade (Rengine = 1.9 um). 


J. Therm. Sci., Vol.25, No.1, 2016 


@ "as casting" blade (Raes = 8.49 um), representing a 
very rough degree of finish (Rdengine= 2.8 um), lower 
than the current standards, but very economic. 

Results obtained with a rough-meter for the three con- 
ditions are reported in Fig. 2. 

The values of dimensionless parameters characterizing 
the different surface finishing are reported in Table 2. Since 
k;` is a local parameter, useful for the description of the 
local influence of the roughness along the pressure and 
the suction side of the blade, the Reynolds number based 


kU, 


on the equivalent sand grain diameter ( Re, = ) is 
here reported to give an overall characterization of the 
different roughness levels. Indeed, for this parameter the 
threshold of Re, = 100 is commonly considered to identi- 


Alpha = 45° Beta = 30° 


OA=NW 3 


O=NWATOANOOSB233 E 


Alpha = 45° Beta = 30° um 


Alpha = 45° Beta = 30° um 


Fig. 2 Roughness distributions for the three conditions (Ra 
increases from top to bottom) 


Table 2 Roughness Reynolds number of the blades tested [A]: 
Adams; [S]: Schäffler 


Re = 70000 Re = 300000 
Rek Rek 
[A] [S] [A] [S] 
“super polished” 3 5 14 21 
“as drawing” 20 30 84 127 
“as casting” 29 44 124 188 
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fies the roughness level influencing the boundary layer 
evolution, Schlichting [24]. 


Results and discussion 


The aerodynamic loading distributions for the high 
and the low Reynolds numbers tested (Re = 300000 and 
Re = 70000) for the steady and the unsteady incoming 
flow conditions are plotted in Fig.3, in the form of the 
static pressure coefficient. For clarity, data reported are 
referred only to the super-polished roughness configura- 
tion. Black and blue curves represent the distributions for 
the steady and the unsteady inflow cases, respectively. 
The blade aerodynamic loadings showed in Fig. 3 for the 
different flow conditions are almost completely supe- 
rimposed, both in the pressure side and in the suction 
side of the blade. The peak velocity on the suction side 
occurs close to x/C, = 0.45, thus a mid-loaded distribu- 
tion characterizes the tested blade. The adverse pressure 
gradient generated by the blade doesn’t give rise to 
boundary layer separation in any flow condition, as re- 
vealed by the absence of changes of slope and plateaus in 
the C, distributions. Indeed, as shown by Lengani and 
Simoni [27] in their measurements on the same test case, 
the high freestream turbulence intensity level keeps the 
boundary layer attached to the wall even in the most un- 
favorable condition of steady inflow at Re = 70000. All 
these considerations referred to the super-polished blade 
can rightly be extended to the as casting and as drawing 
blades, since loading distributions measured in these 
cases well fit those shown in Fig.3. Indeed aerodynamic 
loading distributions are not influenced by the surface 
roughness of the blade. 

The total pressure loss coefficient distributions down- 
stream the "super polished", the "as casting" and the "as 
drawing" blades for the steady incoming flow condition 


1.6 


1.2 


—€— Re=300000 Steady 
—©6— Re=70000 Steady 

—¢— Re=300000 Unsteady 
—©— Re=70000 Unsteady 


0.4 


Fig. 3 Blade aerodynamic loadings 
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are plotted in Fig.4 for Re = 70000 (above) and Re = 
300000 (below) conditions. The pressure side is located 
on the left while the suction side is on the right. The 
overall profile losses coefficient œw of the three roughness 
configurations tested are instead summarized in Table 3. 
Data are made dimensionless by means of the total pres- 
sure loss coefficient measured at Re = 300000 for the 
super polished Ra level. 

For the lower Reynolds number condition Re=70000 
no significant variations can be observed in the distribu- 
tions of C, for the three roughness configurations tested. 
The small deviations between the different cases meas- 
ured experimentally fall amply within the uncertainty 
range of the pressure transducers. Results reported in Fig. 
4 suggest that, for the present steady inflow condition 
characterized by Re = 70000, the profile losses are not 


Re=70000 
0.2 T T T I T T T I T 


— Ra=8.49 um 
— Ra=5.70 um 
—— Ra=0.95 um 


0 0.2 0.4 0.6 0.8 1 
ylg 


Fig. 4 Total pressure losses coefficient distributions in the 


wake region of the three blades tested for steady in- 
flow at Re=70000 (above) and Re=300000 (below). 


Table 3 Overall profile loss coefficients for steady inflow 


Re=70000 Re=300000 
0/Wret O/ ret 
“super polished” 1.705 1.000 
“as drawing” 1.695 1.009 
“as casting” 1.700 1.038 
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influenced by the surface roughness. It is worth noting 
that Re, values reported in Tab. 2 for this flow condition 
are well below the critical threshold of Re, = 100 for all 
the different Ra levels. This means that the roughness 
asperities likely remain confined in the laminar sublayer 
of the boundary layer, and therefore they don’t affect the 
transition process and consequently the profile losses, as 
it will be better documented in the following by means of 
CFD results. On the contrary, for the higher Reynolds 
number (Re=300000), the growth of the total pressure 
defect measured in the blade wake is not negligible as the 
surface roughness increases. In particular, both a growth 
of the height of the peak and a slight increase in the val- 
ues of C, relative to the suction side of the blade can be 
observed in this flow condition. Contrarily, the pressure 
side region of the wake is not influenced. The difference 
between the overall loss values measured for the "super 
polished" and "as casting" blades is greater than the un- 
certainty range. In contrast, the difference of profile 
losses between the "super polished" and the "as drawing" 
cases is contained within the uncertainty range, even 
though varies with a reasonable trend. The similar loss 
values found for the "as drawing" and the "super po- 
lished" blades suggests that both these configurations are 
in the "hydraulically smooth" regime, as also confirmed 
by the Re, values shown in Tab. 2. Conversely, sensibly 
higher values characterize the “as casting” condition. In 
order to assess whether the effects of surface roughness 
on profile losses are altered by the upstream rotor wakes, 
tests have been also carried out in the unsteady flow con- 
dition. The total pressure loss coefficient distributions 
measured downstream the "super polished", the "as cast- 
ing" and the "as drawing" blades are plotted in Fig.5. The 
case Re = 70000 and Re = 300000 are reported on the top 
and bottom side, respectively. The overall dimensionless 
loss coefficient w relative to the three roughness confi- 
gurations tested for the unsteady flow condition are 
summarized in Table 4. Even in this unsteady flow case, 
no significant variations in the distribution of C,, between 
the three roughness configurations have been found in 
the low Reynolds number condition. Consequently, simi- 
lar profile loss coefficients œ characterize the three Ra 
levels, with differences contained within the uncertainty 
range. Thus, also in the presence of the incoming wakes, 
the boundary layer evolution and losses do not seem to 
be influenced by the surface roughness at Re = 70000. 
This suggests that the threshold Re, = 100 can be suc- 
cessfully adopted also in the unsteady flow case. Al- 
though the boundary layer transition process and hence 
its conditions is altered by the disturbances introduced by 
the incoming wakes, the roughness asperities continue to 
remain confined in the laminar sublayer without affecting 
profile losses. 
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As already seen in the steady case, for higher Rey- 
nolds number the total pressure defect measured in the 
blade wake grows appreciably as the surface roughness 
increases. This growing affects both the peak and the 
suction side portion of the profile wake. The profile loss 
coefficient œ grows as the surface roughness increases. 
Similar observations raised during the analysis of the 
steady inflow case can be applied here for the unsteady 
condition, further supporting the applicability of the 
threshold level of Re, = 100 to identify the critical Ra 
levels influencing the blade performance also in the un- 
steady case. However, the effects of roughness seem to 
be more marked in the unsteady flow condition. In the 
unsteady flow cases the difference of losses between the 
extreme configurations is nearly doubled with respect to 
the steady inflow condition, but the relative variations of 


Table 4 Overall profile loss coefficients for unsteady inflow 


Re=70000 Re=300000 
0/Oref 0/ ret 
“super polished” 2.395 2.119 
“as drawing” 2.405 2.171 
“as casting” 2.400 2.186 


Re=70000 


—— Ra=8.49 um 

0.16 — Ra=5.70 pm 

— Ra=0.95 um 
0.12 

Š 
0.08 
0.04 
Ss 
0 
0 0.2 0.4 0.6 0.8 1 


0 0.2 0.4 0.6 0.8 1 
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Fig. 5 Total pressure losses coefficient distributions in the 
wake region of the three blades tested for unsteady in- 
flow at Re=70000 (above) and Re=300000 (below). 
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profile losses are quite similar in the two cases. In the 
unsteady case the effects of incoming wakes further con- 
tribute to shift upstream the boundary layer transition 
process, thus making the surface asperities more pro- 
nounced with respect to the local boundary layer thick- 
ness. 


Numerical results 


On the basis of the previous discussion, the most crit- 
ical conditions for the profile losses increase with respect 
to the “super polished” blade surface, corresponds to the 
one obtained with steady state simulations at Re=300000. 
For this reason, the discussion of numerical results pre- 
sented in the following is limited to this flow condition. 
For the numerical analyses, the conversion factor be- 
tween average Ra and equivalent sand grain roughness ks 
was defined, as suggested by Adams et al. [26], by 
k,=5.863 Ra. 

Figure 6 shows the pressure coefficient distribution 
along the airfoil for various roughness values. The 
agreement between calculations and experiments is seen 
to be good for both the suction and the pressure sides. It 
must be noticed that the impact of roughness level on the 
blade loading distributions is almost negligible. A com- 
parison between blade load distributions obtained with 
transitional and fully turbulent calculations (dashed line 
in Fig. 6) highlights the differences. Departures from the 
fully turbulent solutions are visible between x/C, = 
0.1-0.5 on the pressure side were the transitional flow 
suggests the presence of a separation bubble. On the suc- 
tion side a slight pressure recovery occurs near x/C,=0.80 
for Ra = 8.49um and x/C\=0.93 for Ra = 5.70 um and 
0.95 um. These are associated with the change in the 
boundary layer blockage due to the transition onset trig- 
gered by different roughness levels. Distributions of the 
k, parameter for the investigated roughness values are 
shown in Fig.7 where the hydraulically-smooth threshold 
corresponding to k, = 5 is indicated. It can be seen how 
the pressure side is expected to have a lower roughness 
effect since the values of the roughness parameter are 
generally lower than those found on the suction surface, 
mainly due to the lower local Reynolds number. The 
sharp rise of k, detected on the suction side is related to 
flow transition onset which moves upstream while in- 
creasing the roughness level. Only with the higher rou- 
ghness level the roughness height in wall units slightly 
exceed the threshold value, thus corresponding impact on 
the losses can be expected. 

Figures 8 a) and b) show the contours of turbulent ki- 
netic energy on a blade-to-blade surface for the interme- 
diate and the maximum value of surface roughness. Dif- 
ferent locations of the transition onset are clearly indi- 
cated by the production of turbulence in the suction side 
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boundary layer. The surface roughness increase moves 
the transition onset upstream, increasing the turbulence in 
the diffusing flow toward the rear part of the suction side 
and producing a thicker wake downstream of the trailing 
edge (Fig.8 b). 

Near the maximum curvature of the pressure surface, a 
relatively high turbulent kinetic energy, is associated with 
the formation of a shear layer due to the presence of a 
recirculating region which is only marginally affected by 
the roughness increase. The low velocity on the pressure 
surface and the continuous flow acceleration keeps the 
boundary layer thickness much smaller with respect to 
the suction side. 

The anticipated transition that occurs for the highest 

level of Ra results in a growth of the region affected by 
turbulent boundary layer. Consequently the profile losses 
increase. Thus, the anticipated transition induced by high 
Ra levels has negative effects on profile losses. 
The detailed boundary layer distributions well explain 
the increase in the values of the C,, as surface roughness 
increases found mainly on the suction side of the blade, 
which also suggests that surface roughness has not 
enough effects on the pressure side boundary layer to 
affect losses. Figure 9 compares the measured and com- 
puted loss coefficient distributions along the pitch in the 
wake region. The effect of roughness increase on the 
wake shape is negligible going from “super polished” to 
“as drawing” surfaces. A departure from this shape is 
detected for the maximum roughness by both the mea- 
surements and the computations. In particular numerical 
results predict a thicker wake on the suction-side, due to 
the earlier roughness-induced transition detected for Ra = 
8.49 um and the consequent larger fraction of the suction 
side affected by turbulent boundary layer (see Fig.8 b). A 
similar behavior was observed experimentally by Zhang 
and Ligrani [28] due to the different growth of boundary 
layers on the suction and pressure sides when increasing 
the amount of surface roughness. 
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Fig. 6 Measured and computed blade load distributions, 


steady state results, Re=300000. 
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Fig. 7 Distributions of the k,’ parameter for various rough- 
ness values, steady state results, Re=300000. 


(a) Ra=5.7 um 
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Fig. 8 Turbulent kinetic energy contours, steady state results, 
Re=300000. 


Figure 10 allows an estimation of the critical rough- 
ness Reynolds number which results between 100 and 
120 for both the measurements and the computations. 
The fact that this value agrees well with the experimental 
findings of Feindt [29] suggests that the correlation used 
for the equivalent sand grain roughness is reasonable for 
the case under investigation . 

In order to separate the losses due to boundary layer 
transition induced by roughness and the losses due to the 
roughness value itself, numerical analyses were carried 
out with the fully turbulent assumption for the three 
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roughness values as well as for the smooth surface condi- 
tion. Results are presented in Figure 11 in terms of nor- 
malized losses as a function of the exit Reynolds number. 
The loss value at Re = 300000, obtained for the smooth 
surface, was used for the normalization. The “as drawing” 
and “super polished” surfaces result in almost the same 
loss levels over the investigated Reynolds number range 
and are practically coincident with the smooth surface 
condition. As far as the “as casting” roughness is con- 
cerned, same loss level was found at the lowest Reynolds 
number corresponding to cruise conditions for the blade. 
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Fig. 9 Measured and computed loss coefficient distributions 
along the pitch at x/C,=0.30 from trailing edge (steady 
state, Re=300000). 


O/Oref 


0 50 100 150 
Re, 


Fig. 10 Normalized total pressure loss coefficient as a func- 
tion of the roughness Reynolds number (steady state, 
transitional flow, Re=300000). 
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Fig. 11 Normalized total pressure loss coefficient as a func- 
tion of the exit Reynolds number (steady state, fully 
turbulent flow). 
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A departure from the hydraulically smooth limit was 
found only at Re = 300000 as expected. As shown by 
Marconcini et al. [30] for a similar cascade tested in the 
same wind tunnel, when the steady state boundary layer 
on the blade surfaces is almost laminar and attached, the 
time averaged total pressure losses obtained with period- 
ic unsteady inflow conditions, turned out to be quite 
close to the ones obtained in steady state conditions with 
a fully turbulent approach. 


Conclusions 


The aerodynamic performances of a LPT mid loaded 
cascade have been experimentally and numerically sur- 
veyed for three different surface roughness configura- 
tions, for different Reynolds numbers (Re = 70000 and 
Re = 300000) under both steady and unsteady flow con- 
ditions. 

At Re = 70000 results shown that profile losses are not 
affected by the considered roughness range both under 
steady and unsteady flow conditions. Conversely, at Re = 
300000 profile losses grow monotonously as surface 
roughness increases, both in the steady and in the un- 
steady flow conditions. The absolute effect of roughness, 
however, seems more marked in the unsteady flow con- 
dition, but relative variations of profile losses due to 
roughness are quite similar in the two cases. CFD results 
clearly show that the higher profile losses characterizing 
the high Ra level at Re = 300000 is due to an anticipation 
of the transition process of the suction side boundary 
layer. In fact, this losses increase mainly affect the suc- 
tion side part of the wake, while no appreciable effects 
have been observed on the pressure side region, where 
the dimensionless roughness parameter k; has been 
found to be well below the critical threshold level. 

However, the reduction of profile losses of the "super 
polished" blade is not so sensitive to justify the economic 
burden on the production process, while "as casting" 
blade losses are unacceptable in view of an advanced 
design of the LPT blades. 


Acknowledgements 


This work is part of a joint research project between 
GE Avio, University of Genova, and University of Flo- 
rence. The authors gratefully acknowledge GE Avio for 
granting permission to publish this paper, Benedetta 
Broccolo, Cristian Lizzer, and Fabrizio Demaria of GE 
Avio for the fruitful discussions on the results. 


References 


[1] J. Nikuradse. Laws of Flow in Rough Pipes. VDI-For- 


[6] 


[11] 


[12] 


[13] 


[14] 


[15] 


Investigations on the Roughness Effects on the Aerodynamic Performances of LPT Blades 41 


chungsheft 361, Series B, Vol. 4, 1933. (NACA TM 1292). 
K. A. Flack and M. P. Schultz. Review of Hydraulic 
Roughness Scales in the Fully Rough Regime. ASME J. 
Fluid Eng., 132(4): 041203, 2010. 

A. L. Braslow. Review of the Effect of Distributed Sur- 
face Roughness on Boundary-Layer Transition. AGARD- 
R-254, 1960. 

H. P. Hodson, R. J. Howell. The Role of Transition in 
High-Lift Low-Pressure Turbines for Aeroengines. Prog. 
Aerosp. Sci., 41(6):419-454, 2005. 

R. Vázquez and D. Torre, The Effect of Surface Rough- 
ness on Efficiency of Low Pressure Turbines. ASME 
Paper No.GT2013-94200, ASME Turbo Expo, June 3-7, 
San Antonio, TX, USA, 2013. 

M. Vera, X. F. Zhang, H. P. Hodson , N. Harvey. Separa- 
tion and Transition Control on an Aft-Loaded Ultra-High- 
Lift LP Turbine Blade at Low Reynolds Numbers: High- 
Speed Validation. ASME J. Turbomach., 129(2): 340-347, 
2007. 

R. J. Boyle and R. G. Senyitko, Measurements and Pre- 
dictions of Surface Roughness Effects on Turbine Vane 
Aerodynamics. ASME Paper No.GT2003-38580, ASME 
Turbo Expo, June 16-19, Atlanta, GA, USA, 2003. 

B. Aupoix, P. R. Spalart. Extensions of the Spalart—All- 
maras Turbulence Model to Account for Wall Roughness. 
Int. J. Heat Fluid Flow, 24: 454-462, 2003. 

B. Aupoix. Roughness Corrections for the k-œ Shear 
Stress Transport Model: Status and Proposal. J. Fluid 
Eng., 137(2): 021202, 2015. 

V. Nagabhushana Rao, R. Jefferson-Loveday, P. G. Tucker, 
S. Lardeau. Large Eddy Simulations in Turbines: Influ- 
ence of Roughness and Free-Stream Turbulence. Flow 
Turbulence Combust., 92: 543-561, 2014. 

J. Bellucci, F. Rubechini, M. Marconcini, A. Arnone, L. 
Arcangeli, N. Maceli, V. Dossena, The Influence of 
Roughness on a High-Pressure Steam Turbine Stage: An 
Experimental and Numerical Study. ASME J. Eng. Gas 
Turb. Power, 137(1): 012602, 2015. 

J. P. Bons, A Review of Surface Roughness Effects in Gas 
Turbines. ASME J. Turbomach., 132(2): 021004, 2010. 

F. Satta, D. Simoni, M. Ubaldi, P. Zunino, F. Bertini. Ex- 
perimental Investigation of Separation and Transition Pro- 
cesses on a High-Lift Low-Pressure Turbine Profile Un- 
der Steady and Unsteady Inflow at Low Reynolds Num- 
ber. J. Therm. Sci., 19(1): 26-33, 2010. 

F. Montomoli, H. P. Hodson, F. Haselbach. Effect of 
Roughness and Unsteadiness on the Performance of a 
New Low Pressure Turbine Blade at Low Reynolds 
Numbers. ASME J. Turbomach., 132(3): 031018, 2010. 
M. Marconcini, F. Rubechini, R. Pacciani, A. Arnone, F. 
Bertini. Redesign of High-Lift Low Pressure Turbine 
Airfoils For Low Speed Testing. ASME J. Turbomach., 


42 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


[23] 


134(5): 051017, 2012. 

F. Satta, D. Simoni, M. Ubaldi, P. Zunino, F. Bertini. 
Loading Distribution Effects on Separated Flow Transi- 
tion of Ultra-High-Lift Turbine Blades. J. Propul. Power, 
30(3): 845: 856, 2014. 

D. Simoni, M. Berrino, M. Ubaldi, P. Zunino, F. Bertini. 
Off-Design Performance of a Highly Loaded LP Turbine 
Cascade Under Steady and Unsteady Incoming Flow 
Conditions. ASME J. Turbomach., 137 (7): 071009, 2015. 
A. Arnone. Viscous Analysis of Three—Dimensional Ro- 
tor Flow Using a Multigrid Method. ASME J. Turbo- 
mach., 116(3): 435-445, 1994. 

A. J. Chorin. A Numerical Method for Solving Incom- 
pressible Viscous Flow Problems. J. Comput. Phys., 2: 12 
—26, 1967. 

D. C. Wilcox. Turbulence Modeling for CFD. 2nd ed., 
DCW Industries, Inc., La Cafiada, ISBN 1-928729-10-X, 
1998. 

R. E. Mayle and A. Schultz. The Path to Predicting By- 
pass Transition. ASME J. Turbomach., 119(3): 405-411, 
1997. 

D. C. Wilcox. Formulation of the k-@ Turbulence Model 
Revisited. AZAA J., 46(11): 2823-2838, 2008. 

R. Pacciani, M. Marconcini, A. Fadai-Ghotbi, S. Lardeau, 
M. A. Leschziner. Calculation of High-Lift Cascades in 
Low Pressure Turbine Conditions Using a Three-Equa- 
tion Model. ASME J. Turbomach., 133(3): 031016, 2011. 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


J. Therm. Sci., Vol.25, No.1, 2016 


H. Schlichting. Boundary-Layer Theory, 7th ed., Mc- 
Graw-Hill Inc., New York, ISBN 0-07-055334-3, 1979. 
A. Schaffler, Experimental and Analytical Investigation 
of the Effects of Reynolds Number and Blade Surface 
Roughness on Multistage Axial Flow Compressors. 
ASME J. Eng. Power, 102(1): 5-13, 1980. 

T. Adams, C. Grant, H. Watson. A Simple Algorithm to 
Relate Measured Surface Roughness to Equivalent Sand- 
Grain Roughness. Int. J. Mech. Eng. and Mechatronics, 
1(1): 008, 2012. 

D. Lengani, D. Simoni. Recognition of Coherent Struc- 
tures in the Boundary Layer of a Low-Pressure-Turbine 
Blade for Different Free-Stream Turbulence Intensity 
Levels. Int. J. Heat Fluid Flow, 54: 1-13, 2015. 

Q. Zhang and P. M. Ligrani. Aerodynamic Losses of a 
Cambered Turbine Vane: Influences of Surface Rough- 
ness and Freestream Turbulence Intensity. ASME J. 
Turbomach., 128(3): 536-546, 2006. 

E. G. Feindt. Untersuchungen über die Abhängigkeit 
des Umshclages laminar-turbulent von der Oberfläche- 
nrauhigkeit und der Druckverteilung, Springer-Verlag, 
Berlin, 1956. 

M. Marconcini, R. Pacciani, A. Arnone, F. Bertini. Low- 
Pressure Turbine Cascade Performance Calculations with 
Incidence Variation and Periodic Unsteady Inflow Condi- 
tions. ASME paper No. GT2015-42276, ASME Turbo 
Expo, June 15-19, Montréal, Canada, 2015. 


